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Three topics, the self-consistent resummation of the per- 
turbative expansion for thermal Yang-Mills (YM) theory, 
nonperturbative analysis of Yang-Mills theories in (24-1) di- 
mensions and modification of the Bose distribution for glu- 
ons, are discussed. These topics are related to the mag- 
netic mass or magnetic screening in the quark-gluon plasma. 
The value of the magnetic mass is argued to be close to 
e'^N/2n = g'^TN/2-K for an SU{N) gauge theory. An analytic 
calculation of the string tension for YM2+1, which agrees with 
numerical simulations to within 3%, is another result. 



I. INTRODUCTION 

In this talk, I shall discuss some ideas on how a mag- 
netic mass can be dynamically generated in a nonabelian 
gauge theory. The ultimate goal, of course, is to under- 
stand the thermal perturbative expansion of Yang-Mills 
(YM) theory, at least to the same extent as the zero- 
temperature theory. Let me start with some simple, well 
known, observations about the magnetic mass. 

Why do we need a magnetic mass ? 

It is well known that the thermal perturbative ex- 
pansion of Yang-Mills theory suffers from infrared di- 
vergences. One set of such divergences are associated 
with the so-called hard thermal loops these are of 
the "electric type" in the sense that they can be cured 
by taking account of the electric screening mass. The 
hard thermal loop effective action, in fact, gives a gauge- 
invariant definition of the electric mass and eventually 
a systematic way of reorganizing perturbation theory to 
include the electric screening. However, the perturbative 
expansion also contains other, "magnetic type" , of diver- 
gences which can be eliminated only if there is screening 
of magnetostatic interactions, or in other words, if there 
is a magnetic mass 

How can magnetic mass arise in Yang-Mills theories ? 

At a very heuristic level, it is easy to see why mag- 
netic screening is to be expected. Nonabelian gauge 
theories have magnetic monopole-like field configura- 
tions which are thermally excited and such a monopole- 
antimonopole plasma can screen magnetic interactions 
just as the charges of the electrical type can screen elec- 
trostatic interactions. A more detailed argument can 
be made as follows. In the imaginary time formalism. 



with Matsubara frequencies cOn = 27rnT, where T is the 
temperature, the gauge fields have a mode expansion as 
Ai{x,x'^) = yli.„(a?) exp(i27rnT2;"). At high temper- 
atures and for modes of wavelength long compared to 
l/T, the modes with nonzero Matsubara frequencies are 
unimportant and the theory reduces to the theory of the 
LUn = mode, viz., a three (Euclidean) dimensional Yang- 
Mills theory (or a (2-|-l)-dimcnsional theory in a Wick 
rotated version). Yang-Mills theories in three or (2+1) 
dimensions are expected to have a mass gap and this is 
effectively the magnetic mass of the (3-|-l)-dimensional 
theory at high temperature This gives the qualita- 
tive origin and a quantitative first approximation to the 
magnetic mass. The w„ 7^ modes can give small cor- 
rections to this mass. 

Many indirect arguments have been proposed for the 
existence of a mass gap in YM3 or YM2+1. For example, 
we can consider an SU (2) gauge theory spontaneously 
broken to U{1). This is effectively compact electrody- 
namics which has a mass gap. This feature could survive 
in the unbroken theory if one has sufficient smoothness as 
the parameters are relaxed towards the unbroken phase. 
The recent conjectured connexion between supergravity 
on anti-de Sitter (AdS) space and gauge theories on the 
boundary of this AdS space gives another way to, at least 
qualitatively, argue that there should be a mass gap ^ . 

Why is magnetic mass important ? 

The individual Feynman diagrams of thermal pertur- 
bation theory are divergent and therefore, for meaningful 
calculations, we need the screening masses. The best def- 
inition of continuum Yang-Mills theory at zero tempera- 
ture is in terms of its perturbative expansion. There is a 
systematic way to obtain finite meaningful results for this 
case. At least, term by term, the expansion is well de- 
fined; questions of convergence, Borel summability, etc. 
are another matter. Ultimately, apart from questions of 
calculability, there is the question of principle: Does ther- 
mal perturbative YM theory exist at least to the same 
extent as the zero temperature limit exists ? This is 
important since it may be the only way to define ther- 
mal YM theory. (A lattice definition is more difficult 
for the thermal case, especially for nonequilibrium situ- 
ations. If we have a well-defined continuum theory, it 
can be extended, with some labor, to the nonequilibrium 
situations.) 

I shall discuss three approaches: the self-consistent 
resummation of perturbation theory, a nonperturbative 
analysis of the mass gap in YM2+1 (which is also very 
interesting in its own right) and the possibility of modify- 
ing the thermal distribution to take account of screening. 
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There are a number of other approaches, including lat- 
tice based ones, which have been tried and which I shall 
not discuss here (See O. Phihpsen's article for a dis- 
cussion partly complementary to mine.) 



II. SELF-CONSISTENT RESUMMATION 

In this approach, we reorganize the perturbation ex- 
pansion by adding and subtracting a mass term to the 
action 1^,^,1). The general strategy is as follows. We 
write the action as 



S = Sym + m^Sn 



(1) 



Sym is the standard Yang-Mills part of the action and 
Sm is a mass term which must respect gauge invariance. 
A is a parameter which is taken to have a loop expansion, 

viz., A = A(i) -I- A(2) H , A(*) being the i-th loop term. 

Calculations will be done in a loop expansion wherein 
Sym + m?Sm is taken to be the zeroth order term. Of 
course suitable gauge-fixing and ghost terms are to be 
added to the above action as usual. The quantum effec- 
tive action or the generator of the 1-particle-irreducible 
graphs is then calculated as 



r = {Sym + m^Srn) + (71 - A(i))5™ 



(2) 



where the ellipsis stands for higher loop terms, terms with 
field dependences other than the structure given by Sm 
and ghost-dependent and gauge-fixing type terms. 71 is 
the coefficient of Sm generated by the one-loop diagrams 
with the rules given by Sym + ^^Sm, similarly 72 cor- 
responds to the two-loop contribution and so on. The 
7i are functions of m. The parameter m is the entire 
dynamically generated magnetic mass by definition, and 
hence, we need the corrections to vanish. Thus we need 
71 - A(i) =0, 72 - A(2) = 0, etc. These conditions de- 
termine A. Finally, we do not want to change the theory, 
only rearrange and resum various terms. Thus A must 
itself be m^. We thus get 



71 + 72 + 



(3) 



This is a 'gap equation' which determines m to the order 
to which the calculation is performed. 

We have to choose a specific mass term to implement 
this procedure. As mentioned in the introduction, the 
LUn = mode is the important one and hence, as a first 
step, we can use a three-dimensional mass term and do 
a three-dimensional calculation. The mass term which 
Alexanian and I have used is given by 



Sm, 



J dn K{An,An) 



(4) 



where An = ^AiUi, An — ^AiUi, Ai being the gauge po- 
tentials, rii is a (complex) three-dimensional null vector 



of the form rii = (~ cos 9 cos ip — i sin ip, — cos 6 sin ip + 
icosifySmO). dfl — sin 6d9dip and denotes integration 
over the angles rii. K{An, An) is related to the Wess- 
Zumino-Witten action as well as the eikonal for a Chern- 
Simons theory and is mathematically very similar to the 
hard thermal loop effective action. I have discussed else- 
where the full expression for K{An, An) and some of its 
nice geometrical properties . Carrying out the integra- 
tion over rii, rii y 



Sm. — 



Afi^k)A-ik) S,,^ 



(27r)3 2 
I A'iA]Air^^V,,u 



(5) 



Fyfe(fc,g,_p)) 
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d{k + q + p) 



Bq^ - [q ■ /s)2 



q - k q ■ {q + k) 
k"^ {q + fc)2 
k-{q + k) 



hi kj 



-{q ^ k) 



{qtqjkk + qkQikj + q^qkh) 



(6) 



The self-consistent procedure then gives, to one-loop 
order, a gap equation with the solution B 



1.2 



'2^ 



(7) 



where = g^T , g being the (3+l)-dimcnsional coupling 
constant. 

The primary advantage of this approach is that it is 
perturbatively implementable. However, some criticisms 
have also been raised about this approach. First of all, 
similar calculations have been done with different mass 
terms. Buchmuller and Philipsen have considered a mass 
term generated by coupling to a Higgs scalar field |^ . (A 
similar mass term was also suggested by Cornwall [[7|.) 
They find a value which is (3/4)-th of our value. Jackiw 
and Pi have used a mass term of the form Fij{D'^)~^Fij 
and find a complex mass (which is somewhat worrisome) 

. More recently, they have used a mass term quadratic 
in the gauge potentials; the gap equation part of the cal- 
culation is then equivalent to a "unitary gauge" limit of 
the BuchmuUer-Philipsen calculation and gives the same 
value for the gap The question then arises: Does 
it matter whether different mass terms give different an- 
swers ? Could we take the point of view that one work- 
able mass term proves the generation of the magnetic 
mass ? After all, the self-consistent procedure mixes up 
different loop orders when solving the gap equation and 
hence there is no reason why different mass terms have to 
agree to the order calculations have been performed. An 
optimist might say that different mass terms could agree 
in a full calculation. Nevertheless, the fact that one can 
get different values for m is worrisome. 
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Another potential difficulty has been pointed out by 
Cornwall |10|. The residue Z at the pole of the cor- 
rected propagator in our approach is about 180. This 
is a big change from Z = 1 at the tree level and sug- 
gests that one may run into difficulties, perhaps devel- 
oping a tachyonic mass, at the two-loop level. And also, 
the two-loop corrections are not parametrically small in 
terms of powers of coupling. (If they are smaller than 
the one-loop contribution, it has to be simply due to nu- 
merical factors.) This makes it difficult to get a system- 
atic expansion and solution for the gap equation, a point 
which has been emphasized by Jackiw and Pi ||] . On the 
positive side though, there is a recent two-loop analysis 
which seems to show that the two-loop contribution is 
only about 15 — 20% of the one-loop contribution ||Tl|| . 

So where do we stand on this problem at this time ? 
I believe there are two distinct questions here. Some of 
the criticisms do apply to the estimation of the value of 
the magnetic mass. To be more confident, one needs to 
have some independent nonperturbative analysis, a ques- 
tion to which I shall return in the next section. However, 
there is also the question of including magnetic screening 
in thermal perturbation theory. To make the expansion 
well defined one needs to introduce an infrared cutoff in 
a systematic gauge-invariant way. The specific numeri- 
cal value is not important to this question of principle; 
nor is it important for the practical calculation of some 
types of processes. It is not consistent in the course of a 
perturbative calculation to introduce an ad hoc infrared 
cutoff whose value is taken from a lattice or some other 
nonperturbative analysis. The self-consistent approach 
can be useful in such a context. 

Finally, there is the question of how much the higher 
Matsubara frequencies as well as how much electrostastic 
screening can modify the result. This has been analyzed 
in references [|2j; see also the article by P. Petreczky in 
these proceedings. 

There can be one other way of introducing an infrared 
cutoff within a perturbative expansion. This would in- 
volve changing the distribution function for the gauge 
bosons, a possibility which needs to be explored to see if 
it can provide a better systematic scheme. I shall return 
to this in the last section. 



III. NONPERTURBATIVE ANALYSIS OF YM2+1 

As I have mentioned in the introduction, the mag- 
netic mass may be considered as the mass gap of YM3 
or YM2+1 in a Wick rotated version. Thus analysis of 
YM2+1 gives another approach to the magnetic mass. Of 
course, the existence of a mass gap in YM2+1 or more 
generally, any kind of nonperturbative analysis of this 
theory, is also interesting in its own right. Here I shall 
report on a nonperturbative analysis of YM2+1 which 
we have carried out over the last three years [^^-^. 
We do a Hamiltonian analysis, starting from first prin- 



ciples; quite simply, what we try to do is to write the 
Schrodinger equation for YM2+1 and solve it. This is 
notoriously difficult in field theory in general, but in our 
case, we can make significant progress by using some re- 
sults from two-dimensional conformal field theory. 

Let us consider an SU (A^)-gauge theory in the Aq = 
gauge. The gauge potential can be written as Ai = 
—if^Af, i = 1,2, where t° are hermitian N x A^-matrices 
which form a basis of the Lie algebra of SU{N) with 
[r,t^] = ip'"'t'', Tr(ri'') = i(5''^ Wavefunctions for 
physical states must be gauge-invariant and thus they 
are functions on the space of all gauge potentials mod- 
ulo gauge transformations, i.e., on the gauge-invariant 
configuration space C. It is therefore necessary to make 
a change of variables to gauge-invariant quantities. We 
are not just interested in a formal change of variables. 
We will need to calculate explicitly the Hamiltonian in 
terms of the new variables, the Jacobian for the change 
of variables, the volume measure for the inner product 
of wavefunctions, etc., and then construct eigenstates of 
the Hamiltonian. Thus we need a parametrization of the 
gauge fields such that these steps can be explicitly carried 
out. We start by combining the spatial coordinates xi , X2 
into the complex combinations z = Xi —1x2, z = Xi+ix2 
with the corresponding components for the potential A = 
A, = + zA2), A = A, = \{A^-iA2) = -(A,)t. 
The parametrization we use is then given by 

A,_ ^ -d,MM-\ Ag = M'f-^diM^ (8) 

Here M, are complex SL{N, C)-matrices. Such a 
parametrization is possible and is standard in many dis- 
cussions of two-dimensional gauge fields. A particu- 
lar advantage of this parametrization is the way gauge 
transformations are realized. A gauge transformation 
A, -> = g'^A^g + g-^d,g, g{x) e SU{N) is ob- 

tained by the transformation M AI^^^ = gM. The 
gauge-invariant degrees of freedom are parametrized by 
the hermitian matrix H = AI^ AI, which may be thought 
of as elements of SL{N, C) / SU (N) . Physical state wave- 
functions are functions of H. 

I shall now go through the main results of our analy- 
sis; details of the calculations may be found in references 
Jl3|-p^ . Some of the technical details are also explained 
in the article by D. Karabali in these proceedings. 

Result 1 : The volume measure on the gauge-invariant 
configuration space (which is the measure of integration 
for the inner product of states) can be exactly calculated. 
It is given by |l^jl^ 

dfi{C) = rf/i(i7)e2^'^(^) (9) 

In this equation, dn{H) is the Haar measure for 
SL{N,C)/SU{N), given explicitly by d^(iJ) 
(det r) [Sip] with H parametrized in terms of a real field 
ip^ix) and H-^SH = Sip" rak{v)tk- S{H) is the Wess- 
Zumino-Witten (WZW) action for H given by 
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S{H)= / TiidHdH-^) 



I 

12 



(10) 



This calculation of the measure of integration for the 
inner products does not suffer from Gribov or gauge- 
fixing ambiguities; for more details, see [p"3|~p^ and also 
the article by Karabali in these proceedings. (Inciden- 
tally, K{An, An) in the mass term used in the last sec- 
tion is the same as S{H) above with the identification 
z = X ■ n, z = X ■ h and with all fields being three- 
dimensional.) 

The inner product for states |1) and |2), represented 
by the wavefunctions '^i and ^'2, is given by 



(1|2) = J dfj.{H)e^^^^"^ ^1^2 



(11) 



Observables are ultimately matrix elements of operators, 
i.e., they are inner products. By (pi]), this gives a corre- 
lator of the two-dimensional WZW model for H, which 
is a conformal field theory. In other words, we have the 
exact mapping 



Matrix elements of 
YM2+1 



Correlators of a 
hcrmitian WZW model 



Result 2 : Normalizable wavefunctions arc functions of 
the current J = {N/n)d^H R-^. 

In other words, the wavefunctions are more restricted 
than being just functions of H\ they can only depend on 
H via the specific combination in J . This result follows 
from considerations of integrable representations in the 
WZW theory via a mapping between unitary and her- 
mitian models ||l^ . The quantity J may be considered 
as the gauge-invariant definition of a gluon. Since wave- 
functions depend only on J by this result, we can write 
the Hamiltonian in terms of J and functional derivatives 
with respect to J . 

This result is also consistent with the fact that the Wil- 
son loop operator, which provides a complete description 
of gauge- invariant observables, can be constructed from 
J alone as 



W{C) = TrP e" 



j:^_,(Adz+Adz) ^ rp^p ^(tt/ca) J 



Result 3 : The YM2+1 Hamiltonian is given by 

n = T + v 

^2 



(12) 



(13a) 



T - 



2 

TO 



r{u) 



SJ'^iu) 



2e^ 



dJa{x)dJa{x) 



mN 

where to — e'^N/2'K and 

N Sab 



(u — vy 



tt{u — v) 



(13c) 



(14) 



Result 4 : The kinetic energy operator T obviously has 
the ground state wavefunction <I>o = 1. With the inner 
product (0), this is normalizable. The term J°'{'u)jji^ 
in T counts the number of currents J and each J con- 
tributes m = e^N/2Ti to the eigenvalue. In particular 



T r 



e'N 



r 



(15) 



This is again an exact result, indicating that we should 
expect a gap of m = e'^N/2Tr for the gluons. 



Result 5 : The vacuum wavefunction ^'o is given by 



^'o = cxp 



-^1 Baix) 



(m + \fr. 



Ba{y) 



3 J — and higher J — terms 



(16) 



$0 = 1 is the ground state wavefunction for T, but 
it is not an eigenstate for Ti, — T + V . By treating V 
perturbatively, we can obtain a 1/e^- or 1/TO-expansion 
of the full vacuum wavefunction. Summing up all terms 
with two powers of J in log^^o, we get the result 5 p5| . 

The first term in ( |l^ ) has the correct (perturbative) 
high momentum limit. Thus although we started with 
the high to or low momentum limit, the result ( p6| ) does 
match onto the perturbative limit. The higher terms are 
also small for the low momentum limit. 

Result 6 : The expectation value of the Wilson loop is 
given, for the fundamental representation, by 



(M^f(C)) = constant exp [— ct^c] 



(17) 



where Ac is the area of the loop C and cr, which may be 
identified as the string tension, is obtained as 



iV2- 1 



Svr 



(18) 



This is a prediction of our analysis starting from first 
principles with no adjustable parameters. Notice that 
the dependence on and N is in agreement with large- 
iV expectations, with a depending only on the combi- 



nation e^N as iV ^ 00. Formula 



gives the val- 



(13b) ues V^/e2 = 0.345, 0.564, 0.772, 0.977 for = 2, 3, 4, 5 
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There are estimates for cr based on Monte Carlo simula- 
tions of lattice gauge theory. The most recent results for 
the gauge groups SU{2), SU{3), SUii) and SU{5) are 



1 fe^k 



21.2 



a/e/ = 0.335, 0.553, 0.758, 0.966 



m4 



We see that our 



result agrees with the lattice result to within ~ 3%. 

Given this result, we might venture to say that we 
have a viable approach to nonperturbative phenomena 
in YM2+1 or YM3. 

A comment on the accuracy of our calculation is in 
order. Why is the result so good ? The string tension 
is determined by large area loops and for these, it is the 
long distance part of the wavefunction which contributes 
significantly. In this limit, the 3 J- and higher terms in 
(O) are small compared to the quadratic term. The lat- 
ter leads to the result 

Result 7 : To lowest order in a power series expansion 
for the current J ~ {N/T:)dipata, the Hamiltonian can be 
simplified as 



+ 0a(a;)(m2- V2)0„(f)] + ... (19) 



where 4)a{k) — \/ N kk / (2TTm)(pa{k) , in momentum space. 

In arriving at this expression we have expanded the 
currents and also absorbed the WZW-action part of the 
measure into the definition of the wavefunctions, i.e., the 
operator dia) acts on 



(20) 



The above equation shows that the propagating parti- 
cles in the perturbative regime, where the power series 
expansion of the current is appropriate, have a mass 
m = e'^N/2'K. This value can therefore be identified as 
the magnetic mass of the gluons as given by this nonper- 
turbative analysis. 

We thus see that the self-consistent method, which 
amounts to selective summation of terms in the per- 
turbative expansion gives a value « 1.2(e2iV/27r) while 
the above analysis, which starts from a high e^- 
approximation and then connects to the perturbative 
limit gives a value for the magnetic mass. The 

two values are quite close to each other and also agree 
with other estimates. I expect that the actual numerical 
value of the magnetic mass should be very close to this 
number. 

Finally, I would like to give a short intuitive argument 
for why there is a mass gap, which may help to under- 
stand the exact calculations reported above. The cru- 
cial ingredient is the measure of integration in the in- 
ner product (prj). Writing Ai?, AB for the root mean 
square fluctuations of the electric field E and the mag- 
netic field B, we have, from the canonical commutation 
rules [Ef,A)] = AE AB ~ fc, where k is the 

momentum variable. This gives an estimate for the en- 
ergy 



2 V AB2 



(21) 



For low lying states, we minimize E with respect to AB^ , 
AB"^-^ ~ e'^k, giving £ ^ k. This corresponds to the 
standard photon. For the nonabelian theory, this is in- 
adequate since {TL) involves the factor e2^-S(«). In fact, 

{n) - / dii{H)e^''^^"^ Ue^E^+B^/e^) (22) 



In terms of B, the WZW action goes like S{H) w 
[-(7V/7r)i / B{l/k'^)B + ...]] we thus see that B follows a 
Gaussian distribution of width AS^ « irk^ /N , for small 
values of k. This Gaussian dominates near small k giv- 
ing AB'^ ~ k^{'K/N). In other words, eventhough £ is 
minimized around AB^ ~ fc, probability is concentrated 
around AB^ ~ k'^{Tr/N). For the expectation value of 
the energy, we then find £ ~ e'^N/2'K + 0[k^). Thus 
the kinetic term in combination with the measure factor 
^2NS(H) gQ^j^j lead to a mass gap of order e^N. The ar- 
gument is not rigorous, but captures the essence of how 
a mass gap arises in our formalism. 



IV. MODIFYING THE BOSE DISTRIBUTION 

I now return to the question of thermal perturbation 
expansion for Yang-Mills theory. The nonperturbative 
analysis I have outlined can give a good understanding 
of the origin of the mass gap and even a value for it. 
However, for thermal perturbative YM theory, we would 
like a scheme for magnetic screening which is perturba- 
tively implementable. The selfconsistent scheme has this 
virtue, but is rather difficult from a calculational point 
of view. There is also the further difficulty of the two- 
and higher loop contributions not being parametrically 
smaller. In this section I propose an alternative method. 
In the previous analyses, we have concentrated on the 
spectrum, arguing for a mass gap. The thermal distribu- 
tion of gluons, which is the other essential ingredient in 
any thermal calculation, was unchanged. It is well known 
that the 'bad' infrared behaviour of the Bose distribution 
for gluons of low momenta is what causes difficulties, a 
situation similar to Bose-Einstein condensation. In such 
a situation there will be correlations in the ground state. 
Further, the theory of the mode with zero Matsubara fre- 
quency is equivalent to YM3 at zero temperature with a 
modified coupling constant = g^T. Clearly the ground 
state of YM3 has strong correlations in it. Thus the dis- 
tribution of gluons will certainly be modified. The pro- 
posal is then to start with a modified distribution which 
incorporates magnetic screening (say depending on a pa- 
rameter to) and then to determine the parameter by free 
energy calculations. This is somewhat like a variational 
way of determining the distribution function. 

The number of ways of distributing n^, particles of mo- 
mentum k in the Gk states of energy ujk = \k\ is given 



by 
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w[{n^}] = n 



(Gfc+nfc-l)! 
(Gfe - l)!nfe! 



(23) 



The maximization of this with equal a priori prob- 
abihties for the states and subject to the condi- 
tion of fixed total energy leads to the standard 
Bose distribution. In other words, we can deter- 
mine the most probable occupation numbers fife by 
d/duk [logty[{nfe}] - PiYjk '^ki^k ~ ^)] = 0' where the 
term /3(^j.nfeWfe — E) enforces the condition of fixed 
total energy E via the Lagrange multiplier (3 which is 
eventually identified as 1/T. A simple way to incorpo- 
rate magnetic screening is to assign an a priori proba- 
bility (1 — ^-^I'm.yik state and maximize W — 
nfe(l-e"'^/")"'=M^[{nfc}]. This leads to the distribution 
function 



1 



a(w) = 



-/^ - 1 -f a{uS) 



(j/r, 



1 



(24) 



fik has good infrared behaviour with — > as w ^ 
0; further, for lo ^ m, is exponentially close to the 
standard Bose distribution. The above distribution can 
be derived from the density matrix 



-a+6 



(25) 



P 



-Hy 



(26) 



where Hym is the Yang-Mills Hamiltonian, Hy m is simi- 
lar but constructed out of the , b operators and A, B can 
be defined as eigenopertors of the corresponding Hamil- 
tonians 



YM 



A] = -ujA, 



[Hym,B]=-loB (27) 



A, B are a generalization of a, b to the interacting case. 
By doing free energy calculations, I expect to determine 
the parameter m variationally. These calculations are 
currently under way. 
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a\a,b\b are the usual creation-annihilation operators. 
a\a refer to the gluons, b\b are similar, but have no [4] 
direct meaning; they are just a mathematical device to 
obtain ( p^ in a simple way. 

At first glance, there would seem to be many ways 
of modifying the a priori probabilities to get magnetic 
screening. But actually, the choices are rather limited. 
The modification ( p4| ) has the properties: 

1) It has infrared cutoff. 

2) Modification is exponentially small in the ultraviolet, 
so that moment calculations are not affected. 

3) It is perturbatively implementable. 

4) It corresponds to a simple density matrix, so that it is 
calculationally simple. 

5) It has a gauge-invariant generalization to an interact- 
ing theory. 

Regarding the last point, a generalization of ( p5| ) is 
given by 



[5] 

[6] 

[7] 

[8] 
[9] 
[10] 

[11] 
[12] 

[13] 

[14] 
[15] 
[16] 

[17] 



R. Pisarski, Physica A158, 246 (1989); Phys.Rev.Lett. 
63, 1129 (1989); E. Braaten and R. Pisarski, Phys.Rev. 
D42, 2156 (1990); Nucl.Phys. B337, 569 (1990); ibid. 
B339, 310 (1990); Phys.Rev. D45, 1827 (1992); J. 
Frenkel and J.C. Taylor, Nucl.Phys. B334, 199 (1990); 
J.C. Taylor and S.M.H. Wong, Nucl.Phys. B346, 115 
(1990); R. Efraty and V.P Nair, Phys.Rev.Lett. 68, 
2891 (1992); Phys.Rev. D47, 5601 (1993); R. Jackiw and 
V.P. Nair, Phys.Rev. D48 (1993) 4991; J.P Blaizot and 
E. lancu, Phys.Rev.Lett. 70, 3376 (1993); Nucl.Phys. 
B417, 608 (1994) 

A.D. Linde, Phys.Lett. B96 (1980) 289; D. Gross, R. 
Pisarski and L. Yaffe, Rev.Mod .Phys. 53 (1981) 43. 



[2: 

[3] J. Maldacena, [hep-th/9711200| S.S 



Gubser, I.R. Kle- 
banov and A.M. Polyakov, Phys.Lett. B428, 105 (1998); 

E. Witten, |hep-th/9802"l5q , 9803131. 

A. Billoire, G. Lazarides and Q. Shafi, Phys. Lett. B103 
(1981) 450; O.K. Kalashnikov, JETP Lett. 39 (1984) 405; 
T.S. Biro and B. Muller, Nucl.Phys. A561 (1993) 477; 
U. Heller, F. Karsch and J. Rank, Phys.Lett. B355, 511 
(1995); F. Karsch et al, Nucl.Phys. B 474, 217 (1996); F. 
Karsch, M. Oevers and P. Petreczky, |hep-lat/9807035 . 
V.P. Nair, Phys. Lett. B352, 117 (1995); G. Alexanian 
and V.P. Nair, Phys.Lett. B352, 435 (1995). 
W. Buchmuller and O. Philipsen, Nucl.Phys. B443, 47 
(1995). 

J.M. Cornwall, Phys.Rev. DIO, 500 (1974); ibid, D26, 
1453 (1982). 

R. Jackiw and S.Y. Pi, Phys.Lett. B368, 131 (1996). 
R. Jackiw and S.Y. Pi, Phys.Lett. B403,297 (1997). 
J.M. Cornw all, Phys.Rev. D5 7, 3694 (1998). 

F. Eberlein, |h"ep-ph/98044"60. 



A. Patkos, P. Petreczky and Zs. Szep, Eur. Phys. C5, 337 
(1998). 

D. Karabah and V.P. Nair, Nucl.Phys. B464 (1996) 135; 
Phys.Lett. B379 (1996) 141; Int. J. Mod. Phys A12 
(1997) 1161. 

D. Karabah, C. Kim and V.P. Nair, Nucl.Phys. B524, 
661 (1998). 

D. Karabah, C. Kim and V.P. Nair, |hep-th/9804132| (to 
be published in Phys. Lett.B). 

K. Gawedzki and A. Kupiainen, Phys.Lett. B215 (1988) 
119; Nucl.Phys. B320 (1989) 649; M. Bos and V.P. Nair, 
Int.J.Mod.Phys. A5 (1990) 959. 

M. Teper, Phys. Lett B311 (1993) 223; O. Philipsen, M. 
Teper and H. Wittig, N ucl.Phys. B469 445 (1996); M. 



Teper, hep-lat/9804008 and references therein 



6 



